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The 2D pair-condensate is characterized by a charge ordered state with a "checkerboard" pattern in 
the planes and with an alternating superstructure along the c-axis. We find that Coulomb energy 
gain occurs along the c-axis, which is proportional to the measured condensation energy (Uo) and to 
T c : E\ u k, 2(£ a b/a<) + l) 2 J/o ~ ksT c and is due to inter-layer charge complementarity (charge asym- 
metry of the boson condensate) where £ (, is the coherence length of the condensate and ao ~ 3.9A 
is the in-plane lattice constant. The static c-axis dielectric constant e c and the coherence length £ a t 
are also calculated for various cuprates and compared with the available experimental data and the 
agreement is excellent. 
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I. INTRODUCTION 

Charge ordering seems to be a general phenomenon 
in various layered metal oxides [1]. Among these ma- 
terials charge ordered state (COS) is also found in su- 
perconducting cuprates, although its role plays in high 
temperature superconductivity (HTSC) is still a matter 
of considerable debate [2]. Local probes such as scan- 
ning tunnelling microscopy (STM) revealed recently ev- 
idence for a static checkerboard charge pattern with a 
real-space modulation periodicity of 4ao in the vortex 
core of Bi2212 [3] which is a provocative evidence for 
pinned charge stripes [4] . Although the precise period of 
the modulation is still controversial, however, it seems to 
be comparable with the lattice constant [4]. 

The checkerboard COS of cuprates attracted recently 
the attention of several theoreticians as well [5]. The 
checkerboard COS (CCOS) can be understood as the in- 
plane alternation of holes and anti-holes (electron-hole 
pairs) in such a way that a Cooper wave-function is lo- 
calized on the coherence area [7]. Therefore the black 
and white "fields" of the CCOS correspond to partial 
charges of ~ ±0.16e (pinned to lattice sites), that is the 
optimal hole content found in various cuprates [9,8]. Us- 
ing this simple static picture the sum of the ~ — 0.16e 
partial charges (anti-holes, e.g. the black fields) on the 
checkerboard gives the 2e charge of the charge carrier 
quasiparticle (charge sum rule for the CCOS) [7]. 

Starting from a Cu-0 bond oriented "checkerboard" 
charge pattern with the observed 4ao periodicity [3] we 
propose a simple phenomcnological model for explain- 
ing the 3D character of HTSC in cuprates supported 
by calculations. Within our model the width of the 
checkerboard coincides with the superconducting coher- 
ence length £ a {,. Furthermore we assume the alternation 



of the "checkerboard" charge pattern along the c-axis 
(that is normal to the planes) which leads to Coulomb 
energy gain. Without this assumption inter-layer (IL) 
Coulomb instabilty occurs in layered cuprates due to the 
enormous IL repulsion of holes. We would like to study 
the magnitude of direct Coulomb interaction between 
charge ordered square superlattice layers as a possible 
source of pairing interaction. Our intention is to un- 
derstand HTSC within the context of an IL Coulomb- 
mediated mechanism. The IL charging energy we wish 
to calculate depends on the IL spacing (d), the IL di- 
electric constant e c , the hole content p and the size of 
the superlattice. Finally we calculate the static c-axis di- 
electric constant e c (that is normal to the a6-plane) and 
the coherence length ^ a b for various cuprates which are 
compared with the available experimental observations. 



II. CHECKERBOARD CHARGE PATTERN AND 
THE SUPERCONDCUTING ORDER 
PARAMETER 

The supercondcuting order parameter which corre- 
sponds to the model with a checkerboard charge mod- 
ulation in the planes takes the form of 



^(x, y) = n^ 2 [cos( — tt) + cos(— n)] 



ao 



ao 



(1) 



where the x and y coordinates are varying in the range 
of x,y = [0; £(!(,]. The factor no is the maximal value of 
the charge density at the lattice site centers. Eq. (1) is 
displayed for the coherence area in FIG 2. The modula- 
tion of the order parameter corresponds to the real-space 
modulation of the hole density in the supercondcuting 



FIG. 1. The alternating "checkerboard" charge ordered 
state of the hole-anti-hole condensate in the 4ao x 4ao charge 
ordered bilayer superlattice model. Each lattice sites (opened 
and filled circles) correspond to a CuOi unit cell. Note the 
charge asymmetry between the adjacent layers. The bilayer 
can accomodate a pair of boson condensate (4e). The in- 
ter-layer charge complementarity of the charge ordered state 
is crucial for getting inter-layer Coulomb energy gain. 



FIG. 2. The 3D view of the "checkerboard" charge pat- 
tern of the boson condensate corresponding to the 4ao charge 
modulation represented by the order parameter given in Eq. 
(1). x and y coordinates are given in A. The wells and peaks 
correspond to holes and anti-holes, respectively. 



(SC) state. This kind of order parameter is given earlier 
by Alexandrov [6]. 

The order parameter must satisfy the charge sum rule 
for the boson condensate indicating the localization of 
the pair condensate within the coherence area, 
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\^{x,y)\ 2 dxdy. 



(2) 



Another restriction on *$>(x, y) is that its integral over a 
unit cell whose area is ~ (ao/2) 2 must correspond to the 
hole content, 



\Qh(ah)\ — / 
Jo 



(a /2) 2 



\^(x,y)\ 2 dxdy w 0.16. 



(3) 



This equation reflects the lattice site centered localiza- 
tion of the holes and anti-holes and later on leads to a 
simple electrostatic model where the charge modulation 
for simplicity is replaced by point charges centered in the 
center of the holes and anti holes. 

IL Coulomb energy gain occurs only in that case when 
holes in one of the layers are in proximity with anti-holes 
in the other layer (FIG 1 , IL electrostatic complementar- 
ity, bilayer 5x5 (4ao x 4ao model). An important feature 
is then that the boson condensate can be described by an 
IL charge asymmetry. Therefore we assume an alternat- 
ing charge pattern along the c-axis. The IL coupling of 
the boson-boson pairs in the bilayer 5x5 model naturally 
suggests the effective mass of charge carriers m* w 4m e , 
as it was found by measurements [11]. 

The 5x5 model can be generalized to represent &NxN 
coherence area where N is the real space periodicity of 
the superlattice. In the N x N superlattice model the 
hole and anti-hole partial charges at each CuOi lattice 
sites are q — ±4e/N 2 [7]. Important to note that the 
charge sume rule holds for the characteristic bilayer with 



= 4e where q f hote is the 



a coherence area gahole 
partial anti-hole charge at the ith anti-hole lattice site. In 
other words a pair of a boson condensate can be localized 
within a characteristic bilayer depicted in FIG 1. 

The alternation of the charge pattern along the c-axis 
is reflected in the order parameter by the layer by layer 
alternation of cos and sin functions. For example, the 
adjacent layer is described by 



^(x,y) 



1/2 



y 

sin( — 7r) + sin( — w)]. 



ao 



ao 



(4) 



Without any experimental evidences we assume in this 
paper that the IL Coulomb energy gain is the main con- 
tribution to the condensation energy of the supercond- 
cuting state. This seems to be rather arbitrary assump- 
tion, however, we see no clear cut evidence for the ab- 
plane contribution to the condensation energy Uq (given 
per CuOi unit cell) as well. In this paper we study 
the case when HTSC is governed purely by the c-axis 
Coulomb energy gain (potential energy driven supercon- 
ductivity) and therefore there is no relation to the kinetic 
energy driven mechanism proposed by P. W. Anderson 
[15]. 

The condensation energy of a bilayer with a coherence 
area in the planes (Uob) is given as follows 



U Qb = 2(n + 1) 



£qfa 
a 



1 2 



E. 



IL,SC 



(5) 



where E^ L ' is the Coulomb energy gain in the SC state. 
Uct is the experimental condensation energy given per 
unit cell. Eq. (5) is generalized for multilayer cuprates 
introducing n. For single layer cuprates n = 0, for bilay- 



ers n=l, etc. The factor 



is the number 



of CuOi lattice sites in the planes (within the coherence 



TABLE I. The calculated coherence length of the pair con- 
densate given in a using the experimental condensation en- 
ergies of various cuprates and Eq. (8) at optimal doping. 

T c (K) k B T c (mcV) Uo (^eV/u.c.) O^M C7M 



Bi2201 


20 


1.6 


10" 


~ 8 




LSCO 


39 


2.5 


21 6 


~ 7 


5-8 c 


T12201 


85 


7 


100 ± 20 d 


~ 5 




Hgl201 


95 


7.8 


80 - 107 e 


~ 5 


5' 


YBCO 


92 


7.5 


110 s 


~ 3 


3-4 fe 


Bi2212 


89 


7.3 


95 s 


~ 3 - 4 


4- 6 1 



ao 3.88A, Uo is the measured condensation energy of various cuprates 
in \i cV per unit cell at optimal doping. a from [24], b Uo f» 2 J/mol 
from [16,17], c from [10], d [25], e Uo f» 12 - 16 m.T/g from [18,26] and 
£ a6 from [19], d t/ ss 11 J/mol from [20], e Uo ~ 10 J/mol from [20], 
^ from [19], 9 from [20], 1 from recent measurements of Wang et al., 
£a6 K3 23l(~ 5 - 6ao) [12], from STM images of rof. [3] £ ab f» 4a , h 
from [10], £ a (, Zc is calculated according to Eq. (8) and is also given in 
Table I and ^ s the measured in-plane coherence length given in 

ao ~ 3.9A. The notations are as follows for the compounds: Bi2201 is 
Bi 2 Sr 2 CuO e+s , LSCO {La-L.s5Sro.15Cv.Oi), T12201 (Tl 2 Ba 2 Cu0 6 ) , 
Hgl201 (HgBa 2 CuO i+ s), YBCO (YBa 2 CuO r ) and Bi2212 is 
Bi 2 Sr 2 CaCu 2 Os + s ■ 

area). Factor 2 is applied in Eq. (5) because we calculate 
the condensation energy of a bilayer. 

The IL Coulomb energy is in the SC state 



N t N 



gIL,SC 



47re E( 



EEE 

n—l m—2 ij 



(n) (m) 
(n,m) 

r • • 



(6) 



where r. 



(n,m) 



(n,m) 



is the inter-point charge distance and 
> <i/x, where d/L is the IL distance (CuOi plane 



to plane, i ^ j). n,m represent layer indexes (n =/= m). 



and q\ m> are the point charges of holes and anti-holes 
centered at CuOi lattice 

sites in the nth and mth layers. First the summation 
goes within the bilayer up to iV 2 (the real space period- 
icity of the CCOS, N w ^ + 1) then the IL Coulomb 

J 1 ao ' 

interaction of the basal bilayer are calculated with other 
layers along the c-axis in both direction (n = 1, 2) [7]. Ni 
is the number of layers along the c-axis. When Ni — ► oo, 
bulk E! L ' SC is calculated. 



III. THE CONDENSATION ENERGY AND T c 



According to the perfect correlation found between Uob 
and T c the following formula can be derived using Eq. (5), 



U, 



0b 



2(n+l) 



jab 

a 



1 2 



1 



k B T c . 



(7) 



Therefore the bilayer condensation energy Uob can di- 
rectly be related to the thermal motion at T c . Although 
the mechanism of thermally induced depairing is not un- 
derstood yet, it might be due to the destruction of the 
lattice-CCOS interactions. In this respect polarons or 
bipolarons might play a role in HTSC [5]. 




2 



100 



20 40 60 80 
critical temperature (K) 
FIG. 3. The critical temperature (K) at optimal doping as 
a function of the bilayer condensation energy (2(n+ 1)N Uo, 
meV) using Eq. (7). The straight line is a linear fit to the 
data. The slope of the linear fit is Uob/T c « fcs which is a 
strong evidence of Eq. (7). 



In order to test the validity of Eq. (7) we estimate the 
coherence length of the pair condensate using Eq. (7) and 
using only experimental data, 



£ab ~ »0 



k B T c 



2(n+l)C/o 



- 1 



(8) 



The results are given in Table I as ( m a o unit) and 
compared with the available measured <^fc P - The agree- 
ment is excellent which strongly suggests that Eq. (8) 
should also work for other cuprates. 

The expression Eq. (7) leads to the very simple formula 
for the critical temperature using Eq. (5) and a simple 
Coulomb expression for the IL coupling energy E I C L,SC 
[7] (T c « k- B 'El L ' sc ) 



The plot of Uob (the bilayer condensation energy) 
against T c is shown on FIG 3 using only experimental 
data. Remarkably the data points of various cuprates 
with a variety of critical temperature fit to a line and its 
slope U~ob/T c is the Boltzmann constant k B - The aver- 
age value we get is k B ~ 1-3 ± 0.2 x 10~ 23 J/K which is 
remarkably close to the value of k B = 1-38 x 10~ 23 J/K. 
In the rest of the paper we will present further evidences 
in order to show that the agreement is might not be ac- 
cidental. 



T C (A^, d, e c ) 



E 



ATTe e c k B — - — „— r 

n—l m—2 tj ' ij 



N t N 2 (n) (m) 
' < _(n,m) 



(9) 



When Ni — > oo, bulk T c is calculated. T c can also be 
calculated for thin films when Ni is finite. e c can also be 
derived 



2 2 Ni N 2 „(n) (m) 

e \ ^ \ ^ \ ^ Hi Hj 

Aire k B T c ^ 2^ 2^ („, m) 

n—l m—2 tj ij 



(10) 



TABLE II. The calculated dielectric constant e c using 
Eq. (10) in various cuprates at the calculated coherence length 



£ a b of the charge 


ordered state g 


;iven in 


Table I. 






d(A) 


Tc(K) a 




e c 


exp 


Bi2201 


12.2 


20 


8 


9.9 


12 a 


LSCO 


6.65 


39 


7 


11.3 


23 ±3, 13.5 6 


Hgl201 


9.5 


95 


5 


6.5 




T12201 


11.6 


85 


5 


13.0 


11.3 C 


YBCO 


8.5 


93 


3 


19.4 


23. 6 d 



where ^ s tnc estimated in-plane coherence length given in ao ~ 

3.9A unit, d is the CuOi plane to plane inter-layer distance in A, T c is 
the experime ntal critical temperature. e c is from Eq. (10). e^ xp are the 
measured values obtained from the following references: a [21], b [23], 
or from reflectivity measurements [25], uj p PS bbcm^ 1 [27], A c PS 3^im 
[26], c from [25], d from reflectivity measurements: u p Pd 60cm _1 [27], 
A c sd 0.9,um [14]. 

where a c-axis average of e c is computed when Ni — > oo. 

The calculation of the c-axis dielectric constants e c 
might provide further evidences for Eq. (7) when com- 
pared with the measured values [13,23]. In Table II 
we have calculated the static dielectric function e c using 
Eq. (10) and compared with the experimental impedance 
measurements [22,23]. e c can also be extracted from 
the c-axis optical measurements using the relation [25] 
e c (w) = e c (oo)— c 2 / (cjpA 2 .), where £oo and u p are the high- 
frequency dielectric constant and the plasma frequency, 
respectively [28]. c and A c are the speed of light and 
the c-axis penetration depth. At zero crossing e c (u>) = 

1 /2 

and Up = c/(A c e c (oo)). Using this relation we predict 
for the single layer Hgl201 the low plasma frequency of 
w p k 8 cm -1 using e c = e c (oo) w 6.5 (Table II) and 
A c « 8 Mm [26]. 

IV. CONCLUSION 

Our primary result is that the boson condensate in the 
superconducting state can be described by an alternating 
"checkerboard" type of charge pattern along the c-axis 
which leads to inter-layer charge complementarity and to 
Coulomb energy gain. Our proposal is that this gain is 
converted to the condensation energy of the supercon- 
ducting state, although the detailed mechanism of this 
process still remains unclear. Within our model the pair- 
ing glue is provided by inter-layer coupling. This phsyical 
picture naturally explains the variation of T c system by 
system in various conditions (external and chemical pres- 
sure, multilayers, heterostructures etc.). Although the al- 
ternation of the charge pattern along the c-axis is not yet 
seen experimentally the model might be useful for further 
studies in the future. The reason for this is simple: with- 
out the assumption of inter-layer charge complementarity 
the superconducting state would " suffer" from enormous 
inter-layer Coulomb instability (repulsion) which is cer- 
tainly not the case. Local probes with sufficient depth 



resolution in thin films might detect the presence of such 
supermodulation of the charge pattern along the c-axis 
in the future if the phenomenon exists in nature. Un- 
fortunatelly the microscopic behaviour of pairing is still 
puzzling: the possible conversion of the c-axis free energy 
gain to pairing needs further understanding. Anyhow the 
modulation of the charge density along the c-axis theo- 
retically provides the possibilty of better understanding 
HTSC. 
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I. INTRODUCTION 

Charge ordering seems to be a general phenomenon 
in various layered metal oxides [1]. Among these ma- 
terials charge ordered state (COS) is also found in su- 
perconducting cuprates, although its role plays in high 
temperature superconductivity (HTSC) is still a matter 
of considerable debate [2]. Local probes such as scan- 
ning tunnelling microscopy (STM) revealed recently ev- 
idence for a static checkerboard charge pattern with a 
real-space modulation periodicity of 4ao in the vortex 
core of Bi2212 [3] which is a provocative evidence for 
pinned charge stripes [4] . Although the precise period of 
the modulation is still controversial, however, it seems to 
be comparable with the lattice constant [4]. 

The checkerboard COS of cuprates attracted recently 
the attention of several theoreticians as well [5]. The 
checkerboard COS (CCOS) can be understood as the in- 
plane alternation of holes and anti-holes (electron-hole 
pairs) in such a way that the Cooper wave-function is 
composed of the anti-holes and is localized within the 
coherence area [?] . Therefore the black and white " fields" 
of the CCOS correspond to partial charges of ~ ±0.16e 
(pinned to lattice sites), that is the optimal hole content 
found in various cuprates [7,8]. Using this simple static 
picture the sum of the ~ — 0.16e partial charges (anti- 
holes, e.g. the black fields) on the checkerboard gives the 
2e charge of the charge carrier quasiparticle (charge sum 
rule for the CCOS). 

Starting from a Cu-0 bond oriented "checkerboard" 
charge pattern with the observed 4ao periodicity [3] we 
propose a simple phenomcnological model for explain- 
ing the 3D character of HTSC in cuprates supported 
by calculations. Within our model the width of the 
checkerboard coincides with the superconducting coher- 



ence length Furthermore we assume the alternation 
of the "checkerboard" charge pattern along the c-axis 
(that is normal to the planes) which leads to Coulomb 
energy gain. Without this assumption inter-layer (IL) 
Coulomb instabilty occurs in layered cuprates due to the 
enormous IL repulsion of holes. We would like to study 
the magnitude of direct Coulomb interaction between 
charge ordered square superlattice layers as a possible 
source of pairing interaction. Our intention is to un- 
derstand HTSC within the context of an IL Coulomb- 
mediated mechanism. The IL charging energy we wish 
to calculate depends on the IL spacing (d) , the IL dielec- 
tric constant e c , the hole content p and the size of the 
superlattice. Finally we calculate the static c-axis dielec- 
tric constant e c and the coherence length £ f, for various 
cuprates which are compared with the experimental ob- 
servations. 



II. CHECKERBOARD CHARGE PATTERN AND 
THE SUPERCONDCUTING ORDER 
PARAMETER 

The supercondcuting order parameter (OP) which cor- 
responds to the model with a checkerboard charge mod- 
ulation in the planes takes the form of 



V(x,y) = n 



2| COs( 7r) + COs( — 7r)]. (1) 

ao ao 



where the x and y coordinates are varying in the range 
of x,y — [0;£ a fc]. For simplicity the distribution of the 
OP is neglected in the 3rd dimension and a nearly per- 
fect 2D character is attributed to the condensate. The 
3D anisotropy of the condensate is negligible in the su- 
perconducting (SC) state which is reflected by the ra- 
tio of the in-plane and out of plane coherence lengths 



FIG. 1. The alternating "checkerboard" charge ordered 
state of the hole-anti-hole condensate in the 4ao x 4ao charge 
ordered bilayer superlattice model. Each lattice sites (opened 
and filled circles) correspond to a CuOi unit cell. Note the 
charge asymmetry between the adjacent layers. The bilayer 
can accomodate a pair of boson condensate (4e). The in- 
ter-layer charge complementarity of the charge ordered state 
is crucial for getting inter-layer Coulomb energy gain. 



FIG. 2. The 3D view of the "checkerboard" charge pat- 
tern of the boson condensate corresponding to the 4ao charge 
modulation represented by the order parameter given in Eq. 
(1). x and y coordinates are given in A. The wells and peaks 
correspond to holes and anti-holes, respectively. 



iab/^c ~ 10 [9]. The factor n is the maximal value of 
the charge density at the lattice site centers. Eq. (1) is 
displayed for the coherence area in FIG 2. The modula- 
tion of the order parameter corresponds to the real-space 
modulation of the hole density in the supercondcuting 
(SC) state. This kind of order parameter is given earlier 
by Alexandrov [6]. 

The order parameter must satisfy the charge sum rule 
for the boson condensate indicating the localization of 
the pair condensate within the coherence area, 
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\^(x,y)\ 2 dxdy. 



(2) 



Another restriction on *&(x,y) is that its integral over a 
unit cell whith the area of ~ (ao/2) 2 must correspond to 
the hole content, 



\Qh(ah)\ 
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(ao/2) 



\^{x,y)\ 2 dxdy w 0.16. 



(3) 



This equation reflects the lattice site centered localization 
of the holes and anti-holes and leads to a simple electro- 
static model where the charge modulation for simplicity 
is replaced by point charges centered in the center of the 
holes and anti holes (see later). 

IL Coulomb energy gain occurs only in that case when 
holes in one of the layers are in proximity with anti-holes 
in the other layer (FIG 1 , IL electrostatic complementar- 
ity, bilayer 5x5 (4ao x 4ao model). An important feature 
is then that the boson condensate can be described by an 
IL charge asymmetry. Therefore we assume an alternat- 
ing charge pattern along the c-axis. The IL coupling of 
the boson-boson pairs in the bilayer 5x5 model naturally 
suggests the effective mass of charge carriers m* w 4m e , 
as it was found by measurements [10]. 

The 5x5 model can be generalized to represent a N x N 
coherence area where N is the real space periodicity of 



the superlattice. In the NxN superlattice model the hole 
and anti-hole partial charges at each CuC>2 lattice sites 
are q = ±4e/N 2 . Important to note that the charge sume 
rule holds for the characteristic bilayer with a coherence 
area qf hole = 4e where q? hole i s the partial anti-hole 
charge at the zth anti-hole lattice site. In other words 
a pair of a boson condensate can be localized within a 
characteristic bilayer depicted in FIG 1. 

The alternation of the charge pattern along the c-axis 
is reflected in the order parameter by the layer by layer 
alternation of cos and sin functions. For example, the 
adjacent layer is described by 

Without any experimental evidences we assume in this 
paper that the IL Coulomb energy gain is the main con- 
tribution to the condensation energy of the supercond- 
cuting state. This seems to be rather arbitrary assump- 
tion, however, we see no clear cut evidence for the ab- 
plane contribution to the condensation energy Uq (given 
per CuC>2 unit cell) as well. In this paper we study 
the case when HTSC is governed purely by the c-axis 
Coulomb energy gain (potential energy driven supercon- 
ductivity) and therefore there is no relation to the kinetic 
energy driven mechanism proposed by several authors 
[14]. 

The condensation energy of a bilayer with a coherence 
area in the planes (£7q&) is given as follows 



U 0b = 2(n+l) 



£qfc 

a 



+ 1 



IL,SC 



(4) 



where E^ L ' SC is the Coulomb energy gain in the SC state. 
Uq is the experimental condensation energy given per 
unit cell. Eq. (4) is generalized for multilayer cuprates 
introducing n. For single layer cuprates n = 0, for bilay- 

1 2 



ers n = 1, etc. The factor 



jab 

a 



+ 1 



is the number 



TABLE I. The calculated coherence length of the pair con- 
densate given in a using the experimental condensation en- 
ergies of various cuprates and Eq. (8) at optimal doping. 

T c (K) k B T c (mcV) Uo (^eV/u.c.) O^M C7M 



Bi2201 


20 


1.6 


10" 


~ 8 




LSCO 


39 


2.5 


21 6 


~ 7 


5-8 c 


T12201 


85 


7 


100 ± 20 d 


~ 5 




Hgl201 


95 


7.8 


80 - 107 e 


~ 5 


5' 


YBCO 


92 


7.5 


110 s 


~ 3 


3-4 fe 


Bi2212 


89 


7.3 


95 s 


~ 3 - 4 


4- 6 1 



ao ~ 3.88A, Uo is the measured condensation energy of various cuprates 
in \i cV per unit cell at optimal doping. a from [23], b Uo ss 2 J/mol 
from [15,16], c from [9], d [24], c Uo ~ 12 - 16 mj/g from [17,25] and 
£ o6 from [18], d Uo » 11 J/mol from [19], e t/ ~ 10 J/mol from [19], 
^ from [18], 9 from [19], 1 from recent measurements of Wang et al., 
£a& Ri 23l(~ 5 - 6ao) [11], from STM images of ref. [3] £ ab f» 4a , h 
from [9], £ ab Ic is calculated according to Eq. (8) and is also given in 
Table I and ^ s the measured in-plane coherence length given in 

ao ~ 3.9A. The notations are as follows for the compounds: Bi2201 is 
Bi 2 Sr 2 CuO e+s , LSCO {Lax.s5Sro.15Cv.Oi), T12201 (Tl 2 Ba 2 0'uO 6 ), 
Hgl201 (HgBa 2 C'u0 4+ s), YBCO (YBa 2 CuO r ) and Bi2212 is 
Bi 2 Sr 2 CaCu 2 Os + s- 



of CuOi lattice sites in the planes (within the coherence 
area) . Factor 2 is applied in Eq. (4) because we calculate 
the condensation energy of a bilayer. 

The IL Coulomb energy is in the SC state 
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where 
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Q = EE 
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r ■ ■ 
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(6) 



is the inter-point charge distance and 
' i ; ' > diL, where diL is the IL distance (Cu02 plane 



where 

(n.m) 
. ' 

to plane, i ^ j). n,m represent layer indexes (n ^ m). 



and q^ m> are the point charges of holes and anti- holes 
centered at C11O2 lattice sites in the nth and mth layers. 
First the summation goes within the bilayer up to A^ 2 
(the real space periodicity of the CCOS, iV w ^ + 1) 
then the IL Coulomb interaction of the basal bilayer are 
calculated with other layers along the c-axis in both di- 
rection (n = 1,2). Ni is the number of layers along the 
c-axis. When Ni — ► 00, bulk E 1 ^^ is calculated. 



» 



III. THE CONDENSATION ENERGY AND T c 

The plot of Uqi, (the bilayer condensation energy) 
against T c is shown in FIG 3 using only experimental 
data. Remarkably the data points of various cuprates 
with a variety of critical temperature fit to a line and its 
slope Uob/T c is the Boltzmann constant ks- 




20 40 60 80 
critical temperature (K) 
FIG. 3. The bilayer condensation energy (Uob, meV) as a 
function of the critical temperature (K) at optimal doping. 
The straight line is a linear fit to the data. The slope of 
the linear fit is Uob/T c w ks which is a strong evidence of 
Eq. (7). The error bars denote standard deviations estimated 
from various measurements of the condenstaion energy. 



The average value we get is k B w 1.3±0.2 x 1CT 23 J/K 
which is remarkably close to the value of ks = 1-38 x 
10 _23 J/.ff. In the rest of the paper we will present fur- 
ther evidences in order to show that the agreement might 
not be accidental. According to the correlation found be- 
tween Uob and T c the following formula can be derived 
using Eq. (4), 



U ob - 2(n + l) 



Cab _|_ ^ 

a 



U « k B T c . 



(7) 



Therefore the bilayer condensation energy Uob can di- 
rectly be related to the thermal motion at T c . Although 
the mechanism of thermally induced depairing is not un- 
derstood yet, it might be due to the destruction of the 
lattice-CCOS interactions above T c . In this respect po- 
larons or bipolarons might play a role in HTSC [5] . 

In order to test the validity of Eq. (7) we estimate 
the coherence length of the pair condensate derived from 
Eq. (7) and using only experimental data, 



£ab 



ao 



kBT c 



2(n+l)U 



(8) 



The results are given in Table I as ^ b llc (in ao unit) and 
compared with the available measured £„£ p - The agree- 
ment is excellent which strongly suggests that Eq. (8) 
should also work for other cuprates. 

The expression Eq. (7) leads to the very simple formula 
for the critical temperature using Eq. (4) and a simple 
Coulomb expression for the IL coupling energy E* L ' SC 



T c {N,d, e c ) 



e 2 Q 



47re e c fc B 



(9) 



When Ni — > 00, bulk T c is calculated. T c can also be 
calculated for thin films when JVj is finite. e c can also be 
derived 



TABLE II. The calculated dielectric constant e c using 
Eq. (10) in various cuprates at the calculated coherence length 



£ a b of the charge 


ordered state g 


;iven in 


Table I. 






d(A) 


Tc(K) a 




e c 


exp 


Bi2201 


12.2 


20 


8 


9.9 


12 a 


LSCO 


6.65 


39 


7 


11.3 


23 ±3, 13.5 6 


Hgl201 


9.5 


95 


5 


6.5 




T12201 


11.6 


85 


5 


13.0 


11.3 C 


YBCO 


8.5 


93 


3 


19.4 


23. 6 d 



where £%% lc is the estimated in-plane coherence length given in ao ps 
3.9A unit, d is the CuOi plane to plane inter-layer distance in A, T c is 
the experime ntal critical temperature. e c is from Eq. (10). e^ xp are the 
measured values obtained from the following references: a [20], b [22], 
or from reflectivity measurements [24], cu p ps dScm^ 1 [26], A c £S 3/^m 
[25], c from [24], d from reflectivity measurements: u p 60cm _1 [26], 
A c 0.9jum [13]. 



e 2 Q 



47re k B T c 



(10) 



where a c-axis average of e c is computed when Ni —> oo. 

The calculation of the c-axis dielectric constants e c 
might provide further evidences for Eq. (7) when com- 
pared with the measured values [12,22]. In Table II 
we have calculated the static dielectric function e c using 
Eq. (10) and compared with the experimental impedance 
measurements [21,22]. e c can also be extracted from 
the c-axis optical measurements using the relation [24] 
e c (u>) = e c (oo)— c 2 /(WpA^), where e m andw p are the high- 
frequency dielectric constant and the plasma frequency, 
respectively [27]. c and A c are the speed of light and 
the c-axis penetration depth. At zero crossing e c (w) = 

1 /2 

and bj p = c/(A c e c (oo)). Using this relation we predict 
for the single layer Hgl201 the low plasma frequency of 



lu p w 8 cm 1 using e c = e c (oo) 
A c « 8 urn [25]. 



6.5 (Table II) and 



IV. CONCLUSION 

Our primary result is that the boson condensate in the 
superconducting state can be described by an alternating 
"checkerboard" type of charge pattern along the c-axis 
which leads to inter-layer charge complementarity and to 
Coulomb energy gain. Our proposal is that this gain is 
converted to the condensation energy of the supercon- 
ducting state, although the detailed mechanism of this 
process still remains unclear. Within our model the pair- 
ing glue is provided by inter-layer coupling. This phsyical 
picture naturally explains the variation of T c system by 
system in various conditions (external and chemical pres- 
sure, multilayers, heterostructures etc.). Although the al- 
ternation of the charge pattern along the c-axis is not yet 
seen experimentally the model might be useful for further 
studies in the future. The reason for this is simple: with- 
out the assumption of inter-layer charge complementarity 
the superconducting state would " suffer" from enormous 



inter-layer Coulomb instability (repulsion) which is cer- 
tainly not the case. Local probes with sufficient depth 
resolution in thin films might detect the presence of such 
supermodulation of the charge pattern along the c-axis 
in the future if the phenomenon exists in nature. Un- 
fortunatelly the microscopic behaviour of pairing is still 
puzzling: the possible conversion of the c-axis free energy 
gain to pairing needs further understanding. Anyhow the 
modulation of the charge density along the c-axis theo- 
retically provides the possibilty of better understanding 
HTSC. 
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